Four-wave mixing, quantum control and compensating losses in doped 

negative- index photonic metamaterials 
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Negative refractive index metamaterials (NIMs) 
present a novel class of materials that promise a revo- 
lutionary breakthrough in electromagnetics (for review, 
see, e.g., [1]). Nonlinear optics in such materials remains 
so far a less developed branch of optics. The possibility 
of nonlinear electromagnetic responses in such materials 
attributed to the asymmetry of the voltage-current char- 
acteristics of their building blocks was predicted in [2,3]. 
Recent experimental demonstrations of the exciting op- 
portunities to craft nonlinear optical materials with char- 
acteristics exceeding those in natural crystals arc re- 
ported in [4]. Unique nonlinear-optical (NLO) propaga- 
tion effects associated with three- wave (x' 2 ') coupling 
in NIMs, as compared with their well known counter- 
parts in natural materials, were revealed in [5-8]. The 
striking changes in the optical bistability in a layered 
structure including a NIM layer were shown in [9] . A re- 
view of the corresponding theoretical approaches is given 
in [10]. The most detrimental obstacle toward applica- 
tions of NIMs is strong absorption that is inherent to 
this class of materials. The possibility to overcome such 
obstacles based on three-wave optical parametric am- 
plification (OPA) in NIMs was shown in [7, 8] . A great 
deal of technical problems must be solved, however, in 
order to match the frequency domains of negative in- 
dex (NI), strong NLO response and the phase-matching 
to realize such feasibility. Herewith, we propose and ex- 
plore an alternative approach associated with resonant 
four-wave mixing (FWM) nonlinearities embedded 
in NIMs and tailored through quantum control. The pos- 
sibility of compensating losses and manipulating trans- 
parency, refractive index and nonlinear response of the 
NIM sample with control laser(s) is shown. 

The basic idea of the proposed approach is as follows. 
A slab of NIM is doped by four-level nonlinear centers 
[Fig. 1(a)] so that the frequency o>4 falls in the NI do- 



main, whereas all the other frequencies are in the the pos- 
itive index domain. Below, we show the feasibility to pro- 
duce the transparency and amplification for the signal 
wave at W4 controlled by two lasers at u>i and 0*3. These 
three fields also generate an idler at o>2 = o>3 + o>i — 0*4 , 
which experiences either ordinary, population-inversion, 
or Raman amplification provided by the driving field at 
lo\ and controlled by another driving field at 0*3. The 
amplified idler contributes back to o>4 = o>3 + u)\ — o>2 
through FWM which leads to strongly enhanced OPA. 
We assume that the wave vectors of all waves, k 7 , arc 
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Fig. 1. Scheme of quantum-controlled FWM interaction 
(a) and coupling geometry (b). L04 is signal frequency, 0*2 
is idler, and o>i and o>3 are control fields. 71(0*4) < 0. 

co-directed, which is required for phase-matching. Since 
only o>4 experiences negative refraction and all other fre- 
quencies are in the positive index domain, energy flow at 
0J4 is counter-directed against other waves. Hence, this 
signal wave must enter the slab from its opposite side at 
z = L, which is the exit facet for all other waves entering 
the slab at z=Q [Fig. 1(b)]. Correspondingly, the signal 
exits the slab at z=0. Such backwardness of the signal 
wave leads to counterintuitive, distributed feedback-like 
behavior of the OPA process in NIMs described in [7, 8] . 
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The equations describing the coupled waves at W4 and 
L02 are as follows 

(3) 

da^/dz — — a| exp[iAfcz] + (0:4/2) <Z4, (1) 
da,2/dz = 04 exp[iA/cz] — (02/2)02. (2) 
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( % /w4W2/v e 4 e 2/M4M2)(47r/c)x 2 i- B i £, 3 are 
NLO coupling coefficients; Cj and fij are the dielectric 
permittivities and magnetic permeabilities (which are 
negative at W4); Afc = ki + k^ — ft2 — £4; and Oj are the ab- 
sorption or amplification coefficients. Photon fluxes are 
given by S^, 4/^2, 4 = (c/87r)|a2,4| 2 . The amplitudes of 
the fundamental (control) waves E\ and E3 are assumed 
constant along the slab. Transmittance (amplification) 
at L04 is given by the factor r\^ a ~ \a^{0) / a^{L)\ 2 , where 
L is the slab thickness. Note that the signs in (1) are 
opposite to those in ordinary media, which is due to the 
backwardness of the signal wave. The solution to a sim- 
ilar set of equations and its analysis are given in [7,8]. 

Calculations of the optical constants for embedded 
NLO centers (driven by the control fields) with account 
for constructive and destructive quantum interference 
are performed following the density-matrix technique de- 
scribed in [11]. In our simulations, we used the following 
representative values for relaxation rates: energy level 



relaxation rates T„ = 20 x 10 6 s 1 , 



10 6 s 1 ; partial transition probabilities j g i 
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r 1 ; T ln = 10 iO a" 1 . We 
assumed that A2 = 756 nm, A4 = 480 nm. The results 
of numerical simulations for the optical coefficients en- 
tering equations (1) and (2) are shown in Fig. 2. Here, 
S^4 = LU4 — oj m i ; other resonance detunings flj are defined 
in a similar way. The coupling Rabi frequencies are in- 
troduced as Gi = Eidi g /2h and G3 = E 3 d nm /2h. The 
quantities X4,2 are effective linear susceptibilities and 040 
and Xio denote their resonant values when all the driv- 
ing fields are turned off. Figures 2 (a,b) depict changes in 
the absorption and amplification coefficients for the sig- 
nal and idler, respectively, produced by the control fields. 
(Fig. 2 (b) shows the expanded interval corresponding to 
nonlinear interference resonances.) Figure 2 (c) displays 
changes in refractive indices, and Fig. 2 (d) shows wave- 
vector mismatch. The arrow in Fig. 2 (d) points out the 
interval of 2/4, roughly between 0.29 and 2.34, with five 
detunings j/4 for which Afc = (under the assumption 
that the host material does not introduce any additional 
phase mismatch). Figures 2 (e) and 2 (f) show narrow 
resonances in NLO coupling coefficients, 742, with the 
width on the order of other interference resonances. Fig- 
ure 2 proves the feasibility of manipulating local optical 
parameters through nonlinear quantum interference in- 
duced by the control fields. For the used optical transi- 
tions rates, the magnitude of G ~ 10 12 s _1 corresponds 
to control field intensities / of 10 to 100 kW/ (0.1mm) 2 . 
The optimization of the output signal at z = is de- 
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Fig. 2. Nonlinear interference resonances induced by the 
control fields with G\ = G3 = 50 GHz and fii = Q3 = 
2.5Fi g . (2/4 = n 4 /T m i). 



tcrmined by the interplay between absorption, idler gain 
and FWM, with the later depending on the wave vector 
mismatch. This is a multi-parameter problem involving 
sharp resonance dependencies. Results of our numerical 
analysis of the steady-state solutions [11] to the den- 
sity matrix equations and Maxwell's equations for the 
slowly- varying amplitudes in (1) and (2) are shown in 
Fig. 3. The transmission of the host slab in the NI fre- 
quency domain has been set as 10%. Unlike conventional 
media, the output signal for the waves coupled in the 
NIM slab through OPA represents a set of distributed 
feedback-type resonances [7,8]. Such resonant behavior 
can be observed function of the intensity of the 

fundamental fields, the product of the slab length and 
the density of NLO centers, and the resonance offsets 
for the signal and fundamental fields. Figure 3(a) shows 
these narrow transmission resonances. Here, we intro- 
duce the scaled product of the slab length and the den- 
sity number of embedded centers, L/L ra , through the 
resonance absorption length, L ra = a . Figure 3(b) 
displays the second peak in Fig. 3(a) with greater de- 
tail. It shows that the transparency window is on the 
scale of the narrowest (Raman, in this case) transition 
half-width, as a function of detuning, and the resonant 



2 



1 

0.5 

0- 
50 

1.5 
I 1 

CO 
CO 

E 



30 1 



(a) 



CO 



I 0.5 




20 



10 



—Signal T, . 
— Idler ti 2 /1 000 / \ \ 












G r GHz 45 _ 02 y 4 -y 41 y 4 -y 41 2 4 5 g 3 ,GHz 

(e) (f) 



Fig. 3. Laser- induced transmission resonances in the 
NI frequency domain, fii = Sl 3 = 2.5 • Tjg. (a)-(d): 
d = G 3 = 50 GHz. (a)-(c), (e)-(f): 2 = 0. (c)-(d): 
Va = J/41 = 2.5266. (c): maximum is at a±oL = L/L ra = 
Li/Lra ee 37.02, T 4 w 1 at L/L ra = 36.52. (d): inten- 
sity distribution inside the slab: signal - solid line, idler 
- dashed line; L/L ra = 36.52. (e) and (f): L = L\. (e): 
G 3 = 50 GHz. (f): Gi = 50 GHz. 



absorption length, as a function of length. Figure 3(c) 
depicts the same dependence at the optimum resonance 
offset. Figure 3(d) shows that the intensity of the sig- 
nal inside the slab may significantly exceed its output 
value at z=0, which depends on the ratio of the OPA 
and absorption rates. Here, 772(2) = \a-2(z) / a^ 2 . Re- 
markably, the transparency and amplification occur in 
the frequency range where magnitudes of (74,2! are sub- 
stantially less than their resonant values. Figures 3(e) 
and 3(f) display similar resonance dependence on the 
strength of the control fields. The amplification in the 
maximums in Fig. 3 reaches many orders of magnitude, 
which indicates the feasibility of oscillations and, hence, 
the generation of counter-propagating left-handed signal 
and right-handed idler photons. It is known that even 
small amplification per unit length may lead to lasing 
provided that there is a high-quality cavity (or feedback 
resonance) which effectively increases the distance over 
which the amplification occurs. 



According to Fig. 3(c), characteristic values of 
oyoL ~ 10 are required to ensure the transparency 
and gain. Assuming 040 ~ 10~ 16 cm 2 for the reso- 
nance absorption cross-section, which is typical for dye 
molecules, and N ~ 10 19 cm -3 for the density of 
molecules, we obtain that ~ 10 3 to 10 4 cm -1 , and 
the required slab thickness in the range of L ~10 to 
100 fim. At these values, the contribution of the non- 
linear centers in the refraction index is estimated as 
An < 0.5(A/47r)a4 ~ 10~ 2 - 10~ 3 , which essentially 
does not change the linear negative refractive index. 

In conclusion, we propose the compensation of losses 
in strongly absorbing NIMs through embedded tailored 
optical nonlinearitics. Such a possibility is shown with 
a realistic numerical model. We have studied the res- 
onant FWM-bascd OPA in such composite metamatc- 
rials with a negative refractive index at the frequency 
of the signal and a positive index for all other cou- 
pled waves. The strong nonlinear optical response of 
the composite is primarily determined by the embed- 
ded four-level nonlinear centers and, hence, can be ad- 
justed independently. In addition, we have shown the 
possibility of quantum control of the local optical pa- 
rameters, which employs constructive and destructive 
quantum interference tailored by two auxiliary control 
fields. Frequency-tunable transparency windows in the 
negative-index frequency domain, cavity-free generation 
of entangled counter-propagating photons, and the fea- 
sibility of quantum switching in NIMs have been shown. 
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